Spin nematics and magnetization plateau transition in anisotropic Kagome magnets 
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We study S = 1 kagome antiferromagnets with isotropic Heisenberg exchange J and strong easy 
axis single-ion anisotropy D. For D ^ J, the low-energy physics can be described by an effective 
S = f/2 XXZ model with antiferromagnetic Jz ^ J and ferromagnetic J± ~ /D. Exploiting this 
connection, we argue that non-trivial ordering into a "spin-nematic" occurs whenever D dominates 
over J, and discuss its experimental signatures. We also study a magnetic field induced transition 
to a magnetization plateau state at magnetization 1/3 which breaks lattice translation symmetry 
due to ordering of the and occupies a lobe in the B/ J^-J^/ J± phase diagram. 
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Introduction: Magnets with geometrical frustration 
(competition between different spin exchange interac- 
tions caused by the lattice geometry) exhibit many inter- 
esting properties including spin-liquid like low temper- 
ature phases and unusual spin correlations^^. Kagome 
lattice (see Fig ^ magnets provide many examples of 
this, and several realizations with spins ranging from 
S = 1/2 to S' = 5/32i4^ have been experimentally 
studied. On the theoretical side, numerical and ana- 
lytical work suggests that S — 1/2 isotropic Heisen- 
berg antiferromagnet on the kagome lattice is in an un- 
usual phase with an anomalously large density of singlet 
excitations^ at B = 0. At finite B, there also exists evi- 
dence for the presence of a robust magnetization plateau 
state with magnetization pinned to 1 /3 of the saturation 
moment* . 

A particularly interesting example of frustrated mag- 
netism is provided by the Ni^+ based effective S = 
1 Kagome magnet Ni3V208, in which isotropic ex- 
change interactions compete along with sizeable single- 
ion anisotropy terms (and weak Dzyloshinski-Moriya in- 
teractions) resulting in a rich phase diagram in the pres- 
ence of a magnetic fieldSiiS. Motivated in part by this, 
we consider a kagome lattice model with nearest neigh- 
bour antiferromagnetic spin exchange interaction (J > 0) 
between spin S = 1 ions in the presence of an easy-axis 
single-ion anisotropy (D > 0) along the z axis with the 
Hamiltonian 

H ^ jJ^S.-S.-DY^iS^r-BY^S^ , (1) 
{ij) I I 

where {if) refer to nearest neighbour links of the two- 
dimensional kagome lattice. 

Unlike the unfrustrated case, the D term may have im- 
portant effects in frustrated systems even if not very big. 
We therefore study the limit of large D/J and show that 
interesting physics emerges: We show that the ground 
state at S = is a quantum spin nematic associated with 
ordering of < {S'^)'^ > without ordering of the spin itself. 



Upon increasing the field, magnetization plateaus appear 
at specific magnetization values. Of particular interest 
is a plateau at magnetization 1/3 which we show breaks 
translational symmetry. The corresponding plateau tran- 
sition has a number of interesting properties which we 
discuss. 

When D/J is large and positive and B < J, each spin 
is predominantly in the ruz = ±1 states, and we can 
describe the low energy physics in terms of an effective 
Hamiltonian for (pseudo-) spin 5=1/2 variables cr^. Ex- 
plicit calculation to second order m D/J yields the fol- 
lowing effective low energy Hamiltonian in this regimeii: 

(u) (ij) i 

Here, the a are the usual pauli spin matrices, and the 
parameters of H^s are given by J^ ~ 4 J -I- J"^ /D and 
Jx_ « / D; thus, for large D/J we have Jz/Ji_ ~ 
A:D/J+l + 0{J/D). Clearly, the ground-state of this 
pseudospin 5=1/2 XXZ model for small Jz/J± (which 
is not directly related to the physics of our original 5=1 
problem) must be a ferromagnet polarized in the xy 
plane. Below we analyze the large Jz/Ji_ regime (ap- 
propriate for the large D physics of the original 5=1 
model) separately for i? = or small, and B ^ J. 

When B — 0, the dominant diagonal interaction Jz 
leads to frustration since it is impossible to have all pairs 
of neighboring spins pointing anti-parallel to each other 
along the z axis on the kagome lattice. The ground state 
then lives entirely in the highly degenerate minimally 
frustrated subspace with precisely one frustrated bond 
(parallel spins) per triangle, and is selected by the spin- 
exchange dynamics (J_l). This physics in the present 
J_L > case can be understood straightforwardly by 
thinking in terms of variational wavefunctions (as was 
done recentljiiSii^ on the triangular lattice): Since the 
spin-exchange > is unfrustrated, a good variational 
wavefunction for the small Jz/J± ferromagnet is simply 
l^'i?) = Hi I erf = -1-1). Furthermore, a natural description 
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for the state at large Jz/ J± can be obtained by project- 
ing I'^f) to the minimally frustrated subspace described 
above. Since this subspace admits considerable fluctu- 
ations in the values of Uz, such a projected wavefunc- 
tion l^'oo) continues to gain 'kinetic energy' from spin- 
exchange processes, while minimizing the diagonal inter- 
action energy by construction. 

Thus, x-y ferromagnetic order persists even in the large 
Jz/ J± limit at B = 0, and this remains valid for small 
B as well. Moreover, cr^ correlators in l^'oo) are simply 
given by the T = correlations of the classical Ising 
model on the Kagome lattice, and their short-ranged 
naturai^ rules out any co-existing cr^ spin density wave 
order. [The same conclusion has been reached recently in 
other waye^^ and confirmed numerically^^.] What does 
this analysis imply for the original S — 1 magnet? As 
the pseudospin operator (7+ ^ {S^)"^, the xy ferromagnet 
of the pseudospin magnet actually corresponds to an xy 
spin nematic state where < (S*^)^ >^ but < S >— 0. 
Thus, we conclude that spin-1 Kagome magnets with 
strong easy axis anisotropy order into such a spin nematic 
phase with < (<S'+)^ for B — and its immediate 
vicinity. 

The presence of this nematic ordering is one of our 
main conclusions. As a state that breaks the global U{1) 
symmetry of spin rotations about the easy axis this ne- 
matic will have a gapless linear dispersing 'spin' wave 
which will lead to a contribution to the low temper- 
ature specific heat. Further this state will have a non- 
zero finite spin susceptibility for fields both parallel and 
perpendicular to the easy axis. Despite these similarities 
with conventional ordered antiferromagnets there will not 
be any magnetic Bragg spots in neutron scattering as the 
spin itself is disordered. 

In passing we note that the same considerations 
on a triangular lattice again predict nematic ordering 
which coexists with spin density wave ordering of the z- 
component of the spin - this follows directly from the 
arguments above and the results of Ref. 12 13- 17. 

Returning to the Kagome lattice, as we turn on a 
magnetic field B, the magnetization will initially rise 
smoothly with field since the nematic persists for small B. 
As the field is increased to B ^ J, there will be plateaus 
where the magnetization is field independent and fixed to 
specific commensurate values of the magnetization. We 
now show that for a range of B away from B = there 
is such a plateau state at magnetization 1/3 where the 
ground state is a lattice-symmetry broken spin-density 
wave (SDW) state (in which the z component of the spins 
order as in Fig^. 

Working again with the effective XX Z pseudospin 
Hamiltonian we begin in the extreme limit ot J±/ Jz 
by writing B in terms of a reduced field b as B — Jzb and 
noting that the z coupling and field terms in Hcff can 
be combined and rewritten as ^J2ti'^t ^ 26)^, where 
the sum is now over all triangles t of the kagome lat- 
tice. The physics in this (classical) limit is now clear: 
For < 6 < 1, the energy is minimized by having two of 




FIG. 1: (color online) (a) Periodic Kagome lattice and the 
honeycomb net whose bonds pass through the kagome sites. 
In the plaquette ordered state, red honeycomb edges have 
no dimer (cr^ — +1), while green hexagons resonate via the 
ring-exchange process (shown in (c)). In the alternate colum- 
nar state at the same wavevector, dimers cover all red edges 
(a^ — —1) but not green ones, (b) Schematic phase diagram, 
showing the scans I and II discussed in text. 



the spins in each triangle pointing up and one pointing 
down, which yields a magnetization equal to 1/3 of the 
saturation magnetization, while for & > I, the ground 
state magnetization is locked to the saturation value by 
having all spins pointing up. Thus, one expects a magne- 
tization plateau at 1/3 of the saturation magnetization 
in the vicinity of 6 = 0.5, where the energy gap to change 
in magnetization is largest. In this (classical) limit, the 
ground state has extensive degeneracy, as may be easily 
seen by noting that the manifold of low-energy configu- 
rations can be mapped to the perfect dimer covers of the 
honeycomb lattice whose edges pass through the kagome 
sites (with each down spin corresponding to a dimer cov- 
ering the corresponding honeycomb edge) . 

Let us now turn on a small J±. Apart from an unim- 
portant constant shift in energy, the leading non-trivial 
effect of this perturbation is easily seen to arise at third 
order in degenerate perturbation theory and correspond 
to a 'ring-exchange' term which allows flippable hexag- 
onal plaquettes to resonate with amplitude t ~ —J\/ Jz 
(Fign(c)). This quantum dimer model on the honey- 
comb lattice is known to be in a crystalline 'plaquette' 
state that breaks the lattice translation symmetries of 
the honeycomb lattice in order to maximize the num- 
ber of independently flippable plaquettes from which the 
system can gain kinetic energy This implies a ground 
state with long-range density wave order of the and 
of the bond energies (yfcrj + h.c. (Fig^. 

Thus, the plateau state is stable for large finite Jz/J± 
and is therefore expected to occupy a lobe in the B/Jz- 
Jz/J± plane (Fig^. Clearly, the tip of this lobe rep- 
resents a special point along the locus of plateau tran- 
sitions as the vicinity of the tip is distinguished by the 
presence of low energy 'particle-hole' symmetry corre- 
sponding to equal energies for 'quasiparticle' and 'quasi- 
hole' excitations (here quasiparticles and quasiholes are 
distinguished by the sign of the magnetization deviation 
from 1/3 that they induce by their presence). Given 
that the plateau state breaks lattice translation symme- 
try, the transition to the ferromagnet (nematic) at the 
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FIG. 2; Numerical evidence for direct second order transition 
at the particle-hole symmetric tip of magnetization plateau 
lobe (scan I). Here = S'-l°{Q,ujn = 0)/4L^/3, while g™, = 
1 — (m'*)/(3(m^)^) is the standard Binder cumulant of the 
spin-density wave order parameter and SNtot = 0.55atot 



tip and away presents interesting possibilities: Conven- 
tional Landau theory would predict either an intermedi- 
ate phase with both orders present or a first order tran- 
sition. However recent workiSi has shown that Landau 
theory itself can fail in closely related bosonic models 
- the result in such cases is expected to be an unusual 
direct second order phase transition. 

Numerical results: The foregoing provides the motiva- 
tion for our numerical study of the 5=1/2 XXZ model 
at large Jz/J± and finite field B < 0.5 Jz- We use the 
well-documented stochastic series expansion (SSE) QMC 
methodSSi to access the phase diagram. (At large values of 
Jz/J±, some modifications developed recentljsiiSi were 
used to improve the algorithmic efficiency). Most of our 
data is on -L X L samples (where L is number of unit cells) 
with periodic boundary conditions and L a multiple of six 
ranging from 18 to 30 at inverse temperatures (3 ranging 
from 5/J± to 15/ J±. We use standard SSE estimatorsSSi 
to calculate the ferromagnetic stiffness ps, the equal time 
(C"" {q,T = 0) = (o'^(q')ct^/ (—'?))) and static correlators 

{S^"' {q,ujn = 0) = JodrC^'^' (q,r)) ofcr^, and the static 
correlator of the 'kinetic energy' Ki — (cr+cr^ +h.c.)i on 
link I {S^'^' {q,ujn = 0) = Jo drC^''' {q,T)) (here a and 

a refer to the 3 basis sites and six bond orientations in 
a unit cell, and all site and bond types shown in Fig ^ 
are assigned the coordinates of site type when defining 
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FIG. 3: a) Histograms of relative and absolute phase of all 
order parameters, b) Vertical scan (II) showing plateau state 
at J^/Jx = 8.4. 



the fourier transform). 

By analyzing the L and (3 dependence of the Bragg 
peaks at ±Q — ±(— 27r/3, 27r/3) (components refer to 
projections along Tq and Ti (Fig^l) seen in the static 
correlation functions of and Ki , we conclude that spa- 
tial order is established at these wavevectors when ferro- 
magnetism is destroyed in the plateau state; the observed 
wavevector Q is the ordering wavevector of the plaquette 
and columnar states of Fig ^ The static structure fac- 
tors near the onset of the plateau state also reveal the 
presence of an interesting 'dipolar' structure somewhat 
analogous to the dipolar part of dimer correlators in the 
classical honeycomb lattice dimer model2^. These seem 
to simply reflect the local magnetization 1/3 constraint 
imposed by the B and Jz terms in this region of param- 
eter spaceS^ and persist across the transition into the 
ordered state. 

To further probe the nature of the ordering, we also 
measure the statistics of the phases 0„q, Oko of nine 
complex order parameters Vhq = '^aiQi^n — 0) and 
i'Ka = Ka{Q,uJn = 0). From Fig 13 (a), we see that 
that all relative phases are essentially pinned by the en- 
ergetics of the plateau state and there is one indepen- 
dent phase degree of freedom which we take to be 9no 
(9no = 0, ±27r/3 correspond to the three equivalent pla- 
quette ordered states while 6*110 — ti", ±7r/3 correspond to 
the alternative 'columnar' states at the same wavevec- 
tor). Clearly, this overall phase is very weakly pinned 
(compared to the relative phases) even at low tempera- 
tures relatively far from the critical region, but the pres- 
ence of distinct Bragg peaks in the kinetic energy correla- 
tor strongly suggest that the ordering is of the plaquette 
type as simple caricatures of the columnar state contain 
no dimer resonances. 

The statistics of these phases can be interpreted in 
terms of a Landau free energy F written in terms of the 
order parameters V'na, ^pKa- F = Fi(| iV'ifaP) + 

F2, where Fi = ri|V'„aP + ?^2|V'KaP + Ml(IV'naP)^ + 

M2(|V'ifctP)^ + . . . is insensitive to lattice geometry and 
phase information and F2 encodes phase information. 
When Fi develops a minimum at non-zero values of 
its argument, F2 determines the details of the order- 
ing pattern. The most important terms in F2 are 
quadratic invariants which encode the local magnetiza- 
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tion 1/3 constraint on each triangle (imposed by Jz and 
B ~ O.SJ^), and the kinetic energy gain from having 
pairs of neighbouring sites exchanging spins {—J±_). The 
former may be written as a{\^n ■ + |*n • ^*\^) — 
b\'^n ■ with a and b both expected to be positive 
on energetic grounds (we use the dot product notation 

V ■ W = Y^^VaWa, and = (-0110, V'nl,'0n2), ^ = 

(l,e2W3 e^'^VS)^ J = (1,1,1)), while the latter is given 
by cY^K • 0p with c expected to be negative on energetic 

grounds (here = {;il)Ko,^Ki,^K2,^K3,4'Ki,4>Kb) 

and e = (l,e4'^*/3,e2'^*/3,e4'^*/3,l,e2"/3)), addi- 
tion, symmetries also permit a quadratic cross-term 
[{<^K ■ 0)(*; ■ 0*)e-2W3 + . Landau theory thus 
predicts that the relative phases are such that the natural 
order parameters k ■ and ■ have maximum mod- 
ulus and relative phase of 27r/3 (assuming the coefficient 
of the cross-term is negative), and clearly this correctly 
reproduces the phase relationships exhibited by the data 
(Fig 121 (a)). Finally, the absolute phase is expected to 
be chosen by cubic terms of the form Re [(5'„ ■ il)'^] , and 
our data suggests that the renormalized value of the cor- 
responding coefficient is extremely small even away from 
the phase boundary. 

We have also studied the nature of the phase bound- 
ary between the plateau state and the ferromagnet by 
performing several scans, of which we show data for two 
here. The first scan (scan I shown in FigCJ at constant 
B/Jz = 0.5635 was chosen as it intersects the phase 
boundary at a particle hole symmetric point which we 
identify as the tip of the plateau lobe. The sharpness of 
the crossings seen in the plots (Fig[5J of the Binder cu- 
mulant gm and of psL for different sizes strongly suggest 
that we have reached the asymptotic low temperature 
regime and provide indications that the transition is a 
direct second order transition at Jz/Ji_ ~ 7.68 ± 0.02 
with z = 1. From a scaling collapse of these crossing 
curves, we estimate l/u ^ 1.45 ± 0.2, while similar anal- 



ysis for the order parameter m? gives 2/3/i/ ?» 0.75 ± 0.1; 
however, we emphasize that these are estimates and a 
more detailed study with much larger sizes is needed to 
definitively rule out a very weak first order jump or small 
coexistence region and obtain precision values of expo- 
nents. The second, vertical scan (scan II) was performed 
at Jz/ J± = 8.4 primarily to confirm the existence of the 
plateau state over an appreciable range of _B, and yields 
a plateau state for 0.49 < B/Jz < 0.62 (Fig|3(b)). Near 
the transition points, no clear evidence of a first order 
jump is seen, nor is there a well-resolved region with co- 
existing order parameters (further details regarding the 
nature of these transitions will be discussed separately^^). 

Summary: To summarize, we have demonstrated that 
S — 1 kagome antiferromagnets with moderately strong 
single-ion anisotropy of the easy-axis type exhibit an 
interesting spin-nematic state at and in the vicinity of 
-B = 0, as well as a lattice-symmetry broken spin-density 
wave magnetization plateau state at 1/3 magnetization 
for B ~ 0.5 J. We have also presented numerical evi- 
dence that the transition between these is of an unusual 
direct second order type at least at the tip of the plateau 
lobe in the B/J-D/J plane. We hope our work will pro- 
vide some impetus to look for these effects in anisotropic 
S — I kagome magnets. 
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